The H = xp model revisited and the Riemann zeros 
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Berry and Keating conjectured that the classical Hamiltonian H = xp is related to the Riemann 
zeros. A regularization of this model yields semiclassical energies that behave, in average, as the 
non trivial zeros of the Riemann zeta function. However, the classical trajectories are not closed, 
rendering the model incomplete. In this paper, we show that the Hamiltonian H = x{p + £p/p) 
contains closed periodic orbits, and that its spectrum coincides with the average Riemann zeros. 
This result is generalized to Dirichlet L-functions using different self-adjoint extensions of H. We 
discuss the relation of our work to Polya's fake zeta function and suggest an experimental realization 
in terms of the Landau model. 



One of the most promising avenues to prove the Rie- 
mann hypothesis (RH) is to find a self-adjoint operator 
H whose spectrum contains the imaginary part of the 
non trivial Riemann zeros [I] |2] • This idea was suggested 
by Polya and Hilbert in the dawn of the XX century and 
still, one hundred years later, it remains unproved, as 
well as the RH itself (see [3] for a recent review on phys- 
ical approaches to the RH). There are significant hints 
of the validity of the Polya-Hilbert conjecture. Two of 
them are: the Montgomery- Odlyzko law which states, 
that the local statistics of the Riemann zeros is given by 
the Gaussian Unitary Ensemble (GUE) of Random Ma- 
trix Theory, and the formal similarities between counting 
formulas of zeros in Number Theory and energy levels 
in Quantum Chaotic systems. In this web of relation- 
ships, Michael Berry suggested the existence of a classi- 
cal Hamiltonian whose quantum version would realize the 
Polya-Hilbert conjecture 01 ■ This conjectured Hamilto- 
nian must satisfy the following conditions: i) be chaotic, 
with isolated periodic orbits related to the prime num- 
bers, ii) break time reversal symmetry, to agree with the 
GUE statistics and iii) be quasi-one dimensional. These 
conditions were derived from a formal analogy between 
the fluctuation part of the Riemann-Mangoldt formula of 
the zeros of the zeta function and the Gutzwiller formula 
for the fluctuation term of the counting of energy levels 
in a chaotic quantum system. 

In 1999 Berry and Keating showed that the classical 
Hamiltonian Hd — xp fuUfills conditions ii) and iii) but 
not condition i) 5 . The failure of i) is dramatic be- 
cause this Hamiltonian is integrable, and therefore not 
chaotic, and moreover the classical trajectories are not 
closed, which leads naturally to a continuum spectrum. 
Indeed, the Hamiltonian Hd = xp can be quantized in 
terms of the self-adjoint operator H = {xp + px)/2, with 
p = —ihd/dx, and its spectrum is given by the real line 
[51 [7]. In order to obtain a discrete spectrum, out of 
the xp model, Berry and Keating imposed the conditions 
> ix and \p\ > £p, where the minimal length £x, and 
minimal momentum ip span the Planck area ixip — Stt/i 
in phase space. Subject to these conditions, a particle 
with energy E > describes a truncated hyperbola in 



phase space, 

x{t)=£xe\ p{t) = ^e"\ <t<TE = \og^. (1) 

The area bounded by this trajectory, and the x = ix 
and p — £p axes, measured in Planck units, gives the 
semiclassical number of states 

where the constant 7/8 comes from a Maslov phase. 
Rather remarkably, this formula coincides with the 
asymptotic behaviour of the average term in the 
Riemann-Mangoldt formula [1], where E/h is interpreted 
as the height of a non trivial zero. Incidentally, Connes 
also studied the xp Hamiltonian imposing the constraints 
\x\ < A, IpI < A, where A is a cutoff 8 . In the limit 
A — > cx), one obtains semiclassically a continuum spec- 
trum, where the smooth Riemann zeros appear as miss- 
ing spectral lines. However, a more appropiate interpre- 
tation of Connes's result is that Riemann's formula gives 
a finite size correction to the energy levels. Connes's 
regularization were later derived from the Landau model 
of a particle moving in 2D under the action of external 
magnetic and electric fields [S]. 

A fundamental problem of the Berry-Keating model is 
that the classical trajectories are not closed. The particle 
starts at the phase space point {£x,E/£x), and stops at 
the point {E/£p,£p) in a time Te (see eq.Q). The xp 
hamiltonian breaks time reversal, so the particle cannot 
return to its initial position along the time reversed path. 
Berry and Keating suggested different ways to close the 
trajectories, such as the identification of x and —x , and 
p and —p, or the use of symmetries, but no definite con- 
clusion was reached, and consequently, the connection of 
([2| with the Riemann formula could not be put on more 
solid grounds. 

The aim of this letter is to show that the closure prob- 
lem can be solved by a modification of the xp model that 
preserves several of its features, but makes it into a con- 
sistent quantum model. First of all, we shall constrain 
the motion of the particle to the half line £x < x < oo 
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FIG. 1: Classical trajectories given in eqs. (|5| (continuous 
line) and ([TJ (dotted line) 



while the momentum is allowed to take any real value. 
The classical Hamiltonian is defined as 



He] = X 




PGR 



(3) 



where £p is a coupling constant with dimensions of mo- 
mentum. If \p\ » £p, the extra term added to the xp 
Hamiltonian is negligible, but it becomes dominant if 
\p\ << £p, forbidding the particle to escape to infinity, 
since that would cost an infinite energy. This result is 
made clear by the solution of the Hamilton equations 




p = p 



given by 



x{t) 



\Po\ 



^{pl + ei)e-^'-ei 



(4) 



(5) 



p{t) = ±^{pl+£l)e-^^~£l. 



A complete cycle of a classical trajectory can be de- 
scribed as follows (see fig I]). The particle starts at the 
point A — {£x,Po) (with po| ^ ^p)- Then, x increases 
and p decreases monotonically reaching the turning point 
B = {xni{E),£p), where Xra{E) — E/2£p is the maximal 
elongation. After that, the particle moves backwards to 
the point C — {£x,£p/po)i which is attained in a time 



Te = cosh 



E 
2h 



log 



E 



(E » h) 



(6) 



where h = £x£p should not still be identified with Planck's 
constant 2-Kh. At the point C, the particle bounces off, 
meaning that its momentum £p/po becomes po, and the 
cycle repeats itself, with Te being the period. The lat- 
ter process preserves the total energy, and it is analogue 
of the change in the momentum, p — ^ —p of a particle 
hitting a wall. The classical energies are bounded from 
below by the condition \E\ > Ef = 2h. The minimum 
energy correspond to the static solutions x = £x and 
p = ±£p. 



An extra condition on the Riemann dynamics is the 
existence of complex periodic orbits (instantons) with pe- 
riods Tinst,m = Trim (with m an integer) [5]. The orbits 
([l]) of the xp model are periodic in imaginary time, but 
with a wrong period 27ri. After a complex time Ai = itt, 
the position and momenta change sign, which led Berry 
and Keating to suggest the aforementioned identification 
between x and —x, and p and — p, which in any case does 
not close the orbits. This problem does not arise for the 
Hamiltonian ([s]) , which contains complex periodic orbits 
with a period tti, as can be seen from eq.([5]). 

The semiclassical number of states is given by the 
phase space area swept by the particle measured in units 
of 27r/i, and it is given by 



N{E) 



E 

2TTh 

E 
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1 - 

27Th h 
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v/1 - i2h/Ey 



(7) 



0{E-'), 



E 
2h 



» 1. 



This formula agrees with eq.([2]) if h = 2TTh, up to the con- 
stant term, which has not been considered in ([7|). Let us 
now proceed to the quantization of the classical Hamil- 
tonian ([3]). We choose the normal ordering prescription. 



H = x^ 




(8) 



where 1 /p is the ID Green function satisfying p p ^ = 
p~^p = 1, and whose matrix elements are 



P 



':9{y-x) 



(9) 



with 9{x) the Heaviside step function. H acts on a wave 
function ip as 



H'4>{x) — —ix^ 



dx 



|a;2i/,(a;)| 



dy 



e{y - x)y^ip{y) 
(10) 



This operator is hermitean, i.e. {'ipi\H'ip2) = (^^'0i|V'2)j 
if both wave functions satisfy the non local boundary 
condition 

1 f°° 
n£ie'^^i£x)+£p dxx^{x) = 0. (11) 



where g [0, 27r). To derive (11), we have assumed 
that tl^ jx) decays asymptotically faster that x~^^'^. Using 
eq.( 10 1, the Schroedinger equation HipE = EipE becomes 



an integro-differential equation which can be converted 
into a second order differential equation and a boundary 
condition. The solution of both equations yields a unique 
square integrable eigenfunction given by 



(12) 
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FIG. 2: Absolute value wave functions ipEix), given in eq.( 12 1 
(continuous line), and a;~2 + Tr (dotted line). 



FIG. 3: From bottom to top: plot of — log |E(i?)| (Riemann 
zeros), average Riemann zeros, — log |H|j(_E)| (eigenenergies 
of i/ for /i = 27rfi,i? = 7r/4, and -log|H*(£;)| (Polya zeros). 
The cusp represents the zeros of the corresponding equations. 



where K,y{x) is the modified K-Bessel function (the nor- 
malization factor is not included). The asymptotic be- 



haviour of ( 12 ) is given by 




(13) 



where Xin{E) is the maximal length of the classical tra- 
jectory. Ifx << Xm the wave function tPe{x) behaves, up 
to oscilations, as the eigenfunction a;~2 + V of the quan- 
tum Hamiltonian x^px^ . However, ipsix) drops expo- 
nentially in the classical forbidden region (see fig. 2|. The 
hermiticity of H, requires the eigenfunctions il2) to sat- 



isfy the boundary condition (11), which in turn provides 



the equation for the eigenenergies, 



1 j_ 

9 ^ 9 f, 



^7 Ki 



(14) 

All the solutions of this equation will be real, if the Hamil- 
tonian H is, not only hermitean, but also self-adjoint. To 
verify this property we use the von Neumann theorem: H 
is a self-adjoint operator if the deficiency indices and 
n_ coincide ^lUl lllj . These indices are the number of lin- 
early independent solutions of the equations Wtj; = ±ii/). 
Then if n = n+ ~ n_ , the operator H admits infinitely 
many self-adjoint extensions parameterized by matrices 
of the unitary group U{n). In our case we find that 
n+ = n_ — 1, therefore the self-adjoint extensions corre- 
spond to a phase, that can be identified with the factor 
e^'^ appearing in equations (11) and (14). This ends the 



proof of the reality of all the eigenenergies En- 

If 7^ TT, all the eigenenergies are non vanishing and 
form time conjugate pairs {En, —En} with their asso- 
ciated eigenfunctions being related by the time reversal 
transformation ij;-E„{x) = ip*^ ix). If = tt, there is a 
unique state of zero energy Eq = 0, and eigenfunction 
t/^Eoix) oc a:~2 e~'p^/'', while the non zero energy states 
form again time conjugate pairs. The ground state en- 
ergies ±Eq depend strongly on and can be lower or 
higher than the classical value Eq^. 

To fix the value of i), corresponding to the average Rie- 



mann zeros, we use the asymptotic behaviour of eq.(14|, 

E 



„ , , /47r;i\-^ _^ fE 
-Si^) - ( ) e " cos I — 



which vanishes at 

E , E ^ 1 

log = n H — , 

2TTh ^2he 27r 2 ' 



(15) 



(16) 



li h = 27rh and i? = 57r/4, one recovers the semiclassical 
estimates for N{E) given in eqs. ([2| and ([t]). In refer- 
ences [31 [12] , it is shown that a better estimate of the 
average position of the Riemann zeros is obtained equat- 
ing N{E) to a half integer n+\, rather than an integer, 
which in view of eq.(16) yields i5 = 7r/4 (see fig|3|. 

A confirmation of these results comes from a compar- 
ison with Polya's work on the Riemann S-function |13j 
(see also [UH), 

which is an entire and even function in t, whose zeros co- 
incides with the non trivial zeros of C(^ +ii). Polya made 
a Fourier expansion of (17) and truncated it, obtaining 



S*(i) = 47r2(Xo + «(2^) + A'o_^(2^)), (18) 

which is called Polya's Jake zeta function, since it shares 
several properties with S(i). First of all, the zeros of 
S*(i) and S(t), agree in average, as can be seen using 
the asymptotic expansion [2]. 



7r*2 ^t-ie 



t , t 

cos I - log 

2 ^27re 



7tt 



(19) 



This expression vanishes when the argument of the co- 
sine is n -|- 1 , which confirms the aforementioned rule for 
the average location of the Riemann zeros, and in turn 
the choice d — tt/A. A more remarkable fact is that 
all the zeros of S*(t) are real, as was proved by Polya 
using a general theorem on entire functions [13] . This 
theorem can also be applied to prove the reality of all 
the zeros of S^(i?), a result that we obtained using the 
self-adjointness of the operator H . 
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The RH is a particular case of the generalized Riemann 
hypothesis (GRH), which asserts that all the non trivial 
zeros of the Dirichlet L{x, s)-functions, associated to the 
Dirichlet character x, lie on the critical line Res = ^. 
These functions are defined by a series and associated 
Euler product (Res > 1) 



Xin) 



71=1 



n r 



- xip)p 



(20) 



and their analytic extension to the complex plane. x("-) 
are multiplicative arithmetic functions, i.e. x(rim) = 
x(n)x('7i), x('^ + 9™) = x("-):X(l) = Ij where q is the 
modulus of X- i-functions associated to primitive char- 
acters satisfy the functional relation [T3] , 



?(s,x) = 



r( )i(s,x)=exe(l-s,x) 

(21) 

where is the parity and is the sign of a Gaussian 
sum, 



i-x(-i) 



,1/2^ 



nj e 1 



2 ' ^ 

^ ri=i 

(22) 

A L-function is even (odd) if — 0(1). The Riemann 
zeta function corresponds to the trivial character x('^) = 



l,Vn, with = 0, e^j, = 1. Equation (21) yields the 
average location of the zeros of L{x, s) 



t qt 
— log 



1 

"+2 



(23) 



which leads us to the followi ng identification of parame- 
ters in the H model (see eq.(16)). 



E 

7i 



2Trh 



i?=-(3-2a^-2eJ. (24) 



The Riemann zeta function corresponds to the case q — 1, 
for which h = 2Trh, = tt/4. The correspondence ( 24 ) 
implies that the constant h is quantized as a function of 
the modulus of the L-functions, attaining the classical 
limit, /i — > 0, when q oo. 

A physical realization of the Hamiltonian Q is sug- 
gested by the work of reference [5], which showed that 
Hc\ = xp emerges as the effective Hamiltonian of an elec- 
tron moving in the x — y plane, subject to the action of 
a uniform magnetic field B, perpendicular to the plane, 
and an electrostatic potential V{x,y) = Voxy. If F = 0, 
the electron occupies the lowest Landau level which is 
completely degenerate. This degeneracy is broken by the 
potential V{x,y), which in perturbation theory becomes 
the ID Hamiltonian HoS = ujq xp, where ujq = VqP /h ( 
t = yJhc/eB is the magnetic length). The latter Hamil- 
tonian is obtained replacing y — )■ Pp/h in V{x,y). Con- 
sider now that the particle moves in the half-plane x > t 



and that the electrostatic potential is 

{2nl/qf 



V{x, y) = Vox h/ + 



y 



(25) 



Then, the effective Hamiltonian, in the lowest Landau 
level, in units of wq, becomes (|3|, with the identifications 
(.X — £,£p = 27Th/qi and h — 2Trh/q. We expect the 
parameter i9 to arise from an electric field applied at the 
boundary a; = ^ of the system. 

In summary, we have reformulated the Berry-Keating 
xp model in terms of the classical Hamiltonian Hc\ = 
x{p+£p/p) defined on the half-line x > £x, which posseses 
closed orbits and whose semiclassical spectrum agrees 
with the average Riemann zeros. The quantization of 
this Hamiltonian, yields a self-adjoint operator and a 
non local boundary condition parameterized by an angle 
•d. The spectrum of H agrees asymptotically with the 
semiclassical result and the eigenenergy equation is sim- 
ilar to Polya's fake zeta function that approximates the 
Riemann's S function. The construction is generalized to 
the Dirichlet L~ functions, supporting the idea that the 
GRH could have a proof based on a common quantum 
mechanical model. To achieve this goal one has of course 
to find the quantum origin of the fluctuations of the Rie- 
mann zeros. This work suggest two possible scenarios. 
One is to discretize the dynamics as in the Arnold's cat 
map. The other is to modify the Hamiltonian in a non 
trivial way. Further research is required to clarify which 
path is the best. 
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